Superconductivity mediated by charge fluctuations in layered molecular crystals 
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There is no consensus about the mechanism of the superconductivity or the pairing symmetry 
for layered molecular crystals. Applying slave-boson theory to an extended Hubbard model we show 
that for the and f}" crystal structures the superconductivity is mediated by charge fluctuations 
and the order parameter has d^^y symmetry. This is in contrast to the k-(BEDT-TTF)2X family, 
for which theoretical calculations give superconductivity mediated by spin fluctuations and with 
d^2_y2 symmetry. This is the simplest model that can describe the competition between metallic, 
superconducting, insulating, and charge ordered phases that occurs in the 9 and /3" materials. We 
predict several materials that should become superconducting under pressure. 
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The issue of the interplay of superconductivity, mag- 
netism, and charge ordering is relevant to a wide range 
of strongly correlated electron materials. Examples in- 
clude the copper-oxide (high-temperature) superconduc- 
tors [0, colossal magnetoresistance materials heavy 



fermion compounds 
molecular crystals 



vanadium oxides , and organic 
In particular for the cuprate 
superconductors there is controversy about the relative 
importance of charge fluctuations (associated with dy- 
namical "stripes") and antiferromagnetic spin fluctua- 
tions (associated with the Mott insulator which occurs 
when there is an average of one electron or hole for every 
lattice site). For some heavy fermion compounds recent 
experiments support the idea that the superconductivity 
is mediated by spin fluctuations ||^. 

The family k-(BEDT-TTF)2 X § of molecular crys- 
tals have similarities to the cuprates Q including the 
proximity of superconductivity to a Mott insulator in 
the phase diagram. Although there is an average of 
half a hole per molecule the necessary condition of one 
hole per lattice site is met because the molecules are 
paired up (dimerized) within the k- type crystal struc- 
ture. Theoretical calculations ||^ suggest that the super- 
conductivity has dx2-y2 symmetry (as in the cuprates) 
and is mediated by antiferromagnetic spin fluctuations. 
However, there is controversy about whether experiments 
support this In this Letter, we show theoretically 

that the organic superconductors listed in Table || are 
quite different from the n- type materials and the su- 
perconductivity is mediated by charge fluctuations and 
has dxy symmetry. Our results may also be relevant to 
the recent discovery that the quasi-one-dimensional vana- 
dium bronzes /3'-Cuo.65V2 05 and /3-Nao.33V2 05 (which 
are close to a charge ordered insulator) becomes super- 
conducting under pressure [Q. Previously, it was sug- 
gested by Scalapino, Loh, and Hirsch that spin fluc- 
tuations could mediate d^2_j^2 pairing and charge fluctu- 
ations could mediate d^y pairing. 

The materials considered here consist of lay- 
ers of donor molecules [e.g., BEDT-TTF=bis- 
(ethylenedithiatetrathiafulvalene)] alternating with insu- 



lating layers of anions [e.g., X=SF5CH2CF2S03]. For the 
6 and /3" crystal structures the donor molecules are not 
dimerized and so non-interacting electron models (band 
structure calculations) predict a metallic state due to a 
band which is one quarter filled with holes. However, 
some of these materials are insulators at low temper- 
atures. Mott insulators (resulting from the Coulomb 
repulsion between electrons on a single site) only occur 
for a half-filled band. However, the localization of charge 
(and associated insulating behavior) could result from 
charge ordering due to the Coulomb repulsion between 
electrons on neighbouring sites. Indeed, such charge 
ordering is observed in some of these materials and is re- 
flected in a disproportion of charge between neighbouring 
donor molecules (see Ref. [Qj for a brief review of how 
this is determined experimentally). Depending on the 
anion, temperature and pressure, the materials can be 
either a charge ordered insulator, a metal, or a supercon- 
ductor. (A schematic phase diagram of the 9 materials 
is shown in Fig. 1 of Ref. The common feature of 
the superconductors listed in Table I is that in the phase 
diagram the superconductivity occurs in close proximity 
to the insulating and/or charge ordered phase. Five of 
the superconductors have the very unusual property that 
as the temperature decreases the resistivity is increasing 
before entering the superconducting phase. 

The simplest possible strongly correlated electron 
model which can describe the competition among the 
above phases is an extended Hubbard model at quarter- 
filling on a square lattice . The Hamiltonian is 



<ij> ,(T i 
<ij> i<7 



(1) 



where the operator creates an electron in site i with 
spin cr. t is the amplitude for electrons to go from one site 
to a nearest-neighbours one, 1^ is a nearest-neighbours 
Coulomb interaction, and U is the electron-electron in- 
teraction at a given site. is the chemical potential. A 
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further simplification on model can be made consid- 
ering the fact that U >> V,t, so we can (lx U = oo. 
We have previously studied charge ordering within this 
model and here we briefly summarize the main re- 
sults. For V » t, the model has an insulating phase 
with checkerboard charge ordering {i.e., the wavevector 
associated with charge modulation is (7r,7r), see Fig. 3 
m Ref. §). The spins in the charge ordered state are 
antiferromagnetically coupled due to a fourth order ring- 
exchange process. For V/t < {V/t)c ~ 0.69 we find a 
metallic phase with homogeneous charge density while for 
V/t > {V/t)c the system becomes charge ordered. Hence, 
we find a quantum phase transition from a metallic phase 
to a charge ordered phase with a quantum critical point 
at {V/t)c. It should be stressed that this charge order- 
ing instability is not associated with nesting of the Fermi 
surface since at quarter-filling the Fermi surface diameter 
is much smaller than the length of the (tt, tt) wave vector 
Ip^ . We now give a brief outline of our new theoretical 
calculations which show how charge fluctuations near this 
quantum critical point produce superconductivity. Our 
main results are summarised in Fig. |l|. 

First, we extend the SU(2) spin symmetry of (|l|) to 
SU(A^) by allowing the index a to run from 1 to N. We 
have used slave-boson theory |13| combined with a 1/iV 
expansion jl^, where N is considered to be large. This 
type of approach has previously been applied to a wide 
range of strongly correlated electron systems including 
the Kondo model [p"3|-|l5|, heavy fermions the Hub- 
bard model |l7| , and models for the cuprate supercon- 



ductors [g^ 19|. It has been successful in describing the 
physics of the Kondo effect and, for the Hubbard 
model, it can describe the Mott-Hubbard metal-insulator 
transition . The mean-field theory in the slave bosons 
corresponds to the Gutzwiller approximation [0, and 
so the expansion provides a systematic method to 
calculate corrections to this approximation. We briefiy 
outline the main steps of the approach; details can be 
found elsewhere | |l^ , [T8| ,[7t . The condition J7 — > oo pre- 
cludes doubly-occupied sites and so it is convenient to in- 
troduce the following representation [|l3| for the electron 
operators: cj^. = f}„bi, where f}^ represents a fermion at 
site i which carries spin cr, and 6^, is a boson associated 
with the electron charge located at site i. We impose 
the constraint that either a fermion or a boson can be 
at each lattice site, fi^rfL + ^t^i — by introduc- 

ing a Lagrange multiplier, Ai at each lattice site. We 
expand to the next-to-leading order corrections in l/N, 
arriving at a Hamiltonian which describes an effective 
coupled fermion-boson problem: H = + + H^^^ . 
Hf is the result of taking the average of the boson fields: 
h — {hi) and A = (A^), and describes fermions moving in a 
renormalized band with energy dispersion given by ek = 
4^Tk + A - + W^, with Tk = 2(cos(fc,) + cos(fc„)), 
being the Fourier transform of the hopping operator in 



units of the nearest-neighbour hopping t. (n) is the av- 
erage occupation number of the electrons at each lattice 
site. The effect of Coulomb interactions is then two-fold: 
(i) renormalization of the free electron energies by a fac- 
tor 5^ = N/2 — (n), (n) upward shift in the position of 
the band which is given by A = X]k /{ek)(iT'k -I- W). For 
N = 2 and a quarter fiUed-band, (n) = ^, and so the ef- 
fective mass of the quasiparticles is enhanced by a factor 
of 7n*/m = 2. 

The bosonic part of the Hamiltonian, , describes 
the dynamical and spatial fiuctuations about the mean- 
field solution. The field A^ describes fluctuations in the 
no-double occupancy constraint at each lattice site. The 
real charge fluctuations are described by bi. These bo- 
son fields propagate according to, C(q, ifn), where q is 
the momentum of the boson and Vn — 2TmT is a Mat- 
subara frequency with T being the temperature and n 
an integer number. Finally, Hf~'', couples the fermions 
and the bosons so that they propagate according to: 
-D(q, = j^{D°{q,i^ri)~^ - 7r(q, z^„))-i where 7r(q, i^„) 
is the self-energy of the bosons. These bosons originate 
from the electron-electron repulsion and produce inter- 
action between the quasi-particles. In particular, they 
can induce Cooper pairing of quasi-particles. This is 
in analogy to the pairing of electrons due to phonons 
in elemental metals. However, the mechanism of the 
pairing in the present case is the charge fluctuations. 
The scattering amplitude in the particle-particle chan- 
nel between a quasiparticle with momentum k and an- 
other with momentum — k which scatter to k' and — k', 
is denoted r(q = k — k'). It can also be understood as 
the effective potential between the quasiparticles form- 
ing the Cooper pairs. A divergence in the effective in- 
teraction occurs at the charge ordering instability . 
We find that as the ratio V/t is increased, the poten- 
tial varies its shape, developing singularities at (±7r, ±tt) 
when V/t — > {V/t)c « 0.69, at zero temperature. This 
signals the onset of checkerboard charge ordering. 

In order to look for superconducting instabilities near 
the charge-ordering instability we compute Fermi surface 
averages JT^ , pl| of the effective potential, r(q), weighted 
with the different cubic harmonics, into which the effec- 
tive potential can be decomposed. These have different 
symmetries and they read 



gs* (k) = cos{k^) + cos(fcy) 
9d^2_^2 (k) = coa{k^) - cos(fc,y) 
5<ix„(k) sm{k^) sm{ky) 
.9px(k) = sin(fc^). 



(2) 



From these Fermi surface averages, we find that there 
is attraction between the quasiparticles forming Cooper 
pairs with d^y symmetry at V/t > 0.4 and zero tempera- 
ture. Cooper pairing with other symmetries are found to 
be repulsive for all values of V/t. This is in contrast to 
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the d^2_y2 symmetry found for k-(BEDT-TTF)2X within 
renormalized spin-fluctuation calculations Q . 

The finite-temperature phase diagram obtained is 
shown in Fig. |l|. The line separating the metal from 
the charge-ordered phase is defined from the divergence 
at q = (tt, tt) of the quasiparticle scattering amplitude, 
r(q). The dashed line in Fig. [I] is the extension of 
this line, but we do not compute the transition from the 
charge-ordered phase to the superconducting phase. As 
the temperature is lowered it is possible to go from the 
charge ordered state directly into the superconducting 
phase. This re-entrant behaviour might explain the 
most unusual property [dp/dT < 0) of five of the mate- 
rials listed in Table |. 

The qualitative features of the phase diagram, includ- 
ing the dxy pairing symmetry, turn out to be insensi- 
tive to the details of the band structure and the type of 
charge ordering. First, we have changed the shape of 
the Fermi surface by introducing a next-nearest neigh- 
bour hopping along one of the diagonals of the square 
lattice, t' Wc find that varying the ratio, t' /t, in 

the range < t'/t < 1 changes the shape of the Fermi 
surface significantly (see Fig. 5 in Ref. 0) but does not 
destroy the dxy pairing instability. Second, introducing 
a next-nearest neighbours repulsion, V' , changes the mo- 
mentum dependence of the scattering amplitude. Fur- 
thermore, for sufficiently large values of the ratio V /V 
the singularities in r(q) shift from (tt, tt) to (0, ±7r) and 
(±7r,0). This is because for V /V > 1, it is energetically 
more favourable to produce charge ordering either along 
the X or y-direction, rather than along both of them. We 
find that there is still a quantum critical point, and super- 
conducting pairing with dxy symmetry persists, although 
the strength of the effective interaction within the Cooper 
pairs decreases somewhat as compared to the square lat- 
tice case. The robustness of the dxy superconducting in- 
stability can be understood from the fact that the charge- 
ordering instability is not associated with Fermi surface 
nesting. Furthermore, it means our results should also 
be applicable to the /3" materials for which the Fermi 
surface is more complicated than for the materials. 

Proposed experimental tests. (i) We identify sev- 
eral materials that might become superconducting un- 
der pressure. Pressure in the 9 materials decreases 
the hopping parameter t, driving the materials into 
the insulating phase Therefore, 9 materials 

such as 6i-(BETS)2Ag(CN)2, 6'-(BETS)4Cu2(Cl)6 and 9- 
(B0)2C1(H20)3, which are all metallic at ambient pres- 
sure, should become superconducting under pressure. 
These materials would then be located in the metallic 
side of our phase diagram (see Fig. |]), and pressure 
would increase V/t driving them into a superconduct- 
ing state with dxy symmetry before becoming insulating, 
(ii) Polarisation-dependent Raman scattering should be 
done in the superconducting state because it can distin- 
guish d^jj^ and da;2_2^2 states (iii)Measurements of the [12] 



nuclear magnetic resonance relaxation rate and Knight 
shift should be done in the metallic phase for the super- 
conductors in Table I. There should be no enhancement 
of the Korringa ratio. This is in contrast to the large en- 
hancements seen in k-(BEDT-TTF)2X superconductors 
which are close to the Mott insulator Q . 

In conclusion, we have identified 11 molecular super- 
conductors which we predict to have pairing of dxy sym- 
metry due to charge fluctuations. This is based on a 
systematic many-body calculation using slave boson the- 
ory for an extended Hubbard model. Materials with the 
9 and /?" crystal structures are at quarter filling and are 
described by quite different physics from j3, k, and A 
structures which are essentially at half filling. 
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FIG. 1. Phase diagram showing competition between 
metallic, superconducting (SC), and charge ordered phases. 
The symmetry of the Cooper pairs in the superconducting 
phase is d^,,,. The phase diagram is for the extended Hub- 
bard model (defined by the Hamiltonian (^) with an average 
of one hole per two molecules (a quarter-filled band) in the 
limit of infinite Coulomb repulsion energy U for two holes on 
the same molecule. This is the simplest model Hamiltonian 
which can describe the organic superconductors listed in Ta- 
ble I. The vertical axis is the ratio of the temperature T to 
the intermolecular hopping integral t. The horizontal axis is 
the ratio of the nearest-neighbour Coulomb repulsion energy 
V between electrons on neighbouring molecules to t. Note 
that the unconventional superconductivity is found near the 
quantum critical point (at {V/t)c) separating the metallic and 
charge ordered phases. 



4 



TABLE I. Organic superconductors described by our theory. All materials (except the last one) have an average of half 
a hole per donor molecule. This corresponds to a quarter-filled band. Tc is the superconducting transition temperature 
at the given pressure. Five of the materials have the unusual property that the resistivity, p, decreases with increasing 
temperature above Tc suggesting a direct transition from an insulating phase into a superconducting phase {dp/dT < 0). There 
is evidence of charge ordering (CO) and/or an insulating phase (with a metal-insulator transition temperature Tmi) in close 
proximity to the superconducting phase. 6l-(BEDT-TTF)2 I3 is close to a charge ordered insulator as when the anion, I3, is 
replaced with RbZn(SCN)4, CsZn(SCN)4, CsZn(SCN)4 or CsCo(SCN)4, the material becomes a charge- ordered insulator [ p5[ . 
(Y=(C204)3-PhCN). 



Material 


Pressure 


Tc(K) 


Ref. 


dp < 


Pressure 


CO 


Tmi{K) Ref. 


6l-(BEDT-TTF)2l3 


1 bar 


3.6 




2f 


1 


No 










/3"-(BEDT-TTF)2 SF5CH2CF2SO3 


1 bar 


5.2 




27 




No 


1 bar 


Yes 




28 


(BEDT-TTF)2Re04 


4 kbar 


2 




i 

31 


1 


No 


1 bar 


Yes 


77 


30 


/3"-(BEDT-TTF)4Pd(CN)4H20 


7.0 kbar 


1.2 


1 




Yes 


1 bar 




70 


31 


/3"-(BEDT-TTF)4Pt(CN)4H20 


6.5 kbar 


2.0 




2! 




Yes 


1 bar 




120 


25 


/3"-(BEDT-TTF)4H30FeY 


1 bar 


8.5 








No 


1 bar 


Yes 




33 


/3"-(BEDT-TTF)4H30CrY 


1 bar 


5.5 




32 


1 


No 


1 bar 


Yes 




33 


6(-(BETS)2(Cl2 TCNQ) 


3.5 kbar 


1.3 


1 




1 


Yes 


8.5 kbar 




22 


34 


/3"-(BEDO-TTF)2Re04 • H2O 


1 bar 


2-3 








Yes 










6»-(DIETS)2Au(CN)4 


10 kbar uniaxial 
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1 




1 


No 


1 bar 




220 


37 


/3"-(BEDT-TTF)3Cl2 • 2H2O 


16 kbar 


2-3 








Yes 


1 bar 


Yes 


150 


39 
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